MADPH-07-1494 

August 2007 



High energy neutrinos from neutralino annihilations in the Sun 

Vernon Barger/'Q Wai-Yee Keung,^'Q Gabe Shaughnessy/ and Adam Tregre^'y 

^Department of Physics, University of Wisconsin, 
1150 University Avenue, Madison, Wisconsin 53706 USA 
^Physics Department, University of Illinois at Chicago, Illinois 60607-7059 USA 

Abstract 

Neutralino annihilations in the Sun to weak boson and top quark pairs lead to high-energy neu- 
trinos that can be detected by the IceCube and KM3 experiments in the search for neutralino 
dark matter. We calculate the neutrino signals from real and virtual WW, ZZ, Zh, and ti pro- 
duction and decays, accounting for the spin-dependences of the matrix elements, which can have 
important influences on the neutrino energy spectra. We take into account neutrino propagation 
including neutrino oscillations, matter-resonance, absorption, and Vr regeneration effects in the 
Sun and evaluate the neutrino flux at the Earth. We concentrate on the compelling Focus Point 
(FP) region of the supergravity model that reproduces the observed dark matter relic density. 
For the FP region, the lightest neutralino has a large bino-higgsino mixture that leads to a high 
neutrino flux and the spin-dependent neutralino capture rate in the Sun is enhanced by 10'^ over 
the spin-independent rate. For the standard estimate of neutralino captures, the muon signal rates 
in IceCube are identifiable over the atmospheric neutrino background for neutralino masses above 
Mz up to 400 GeV. 



* Electronic address: barger@pheno.physics.wisc.edu 



'''Electronic address: keung @uic.edu| 
Electronic address: gshau@hep.wisc.edu| 
^Electronic address: |tregre@hep.wisc.edu| 



1 



I. INTRODUCTION 



A stable Weakly Interacting Massive Particle (WIMP) of mass of order 100 GeV that 
was produced in the early Universe, thermalized, and froze-out due to the Hubble expan- 
sion provides a natural explanation for observed density of dark matter today [1]. A well 
motivated dark matter (DM) candidate is the Lightest Stable Particle (LSP) of Supersym- 
metry (SUSY) with R-parity conservation . The LSP is nominally the lightest neu- 

tralino (denoted by x?)? a neutral spin-1/2 particle that is a linear combination of gauginos 
(spin-1/2 SUSY companions of the Spin-1 Standard Model bino and Wino) and higgsinos 
(spin-1/2 companions of two spin-0 Higgs bosons) The minimal supergravity model 
(mSUGRA) has minimal SUSY particle content with the exact supersymmetry broken 
gravity jol . With gauge and Yukawa coupling at the Grand Unified Scale (GUT) 7 
the mSUGRA model predictions are given in terms of a small number of parameters 10 
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The consequences of mSUGRA for neutralino dark matter detection have been the subject 



of numerous phenomenological studies that have considered t 



le predictions for forthcoming 



experiments; some recent surveys are given in Refs. [12|, [131, |lj, |15|, [161, |17|, [181, [19]. There 
are three complementary experimental approaches: direct detection via nuclear recoils from 
WIMP scattering, indirect detection via astrophysics experiments wherein WIMP annihila- 
tions give neutrino, gamma ray, positron, and antideuteron signals, and collider experiments 
where the supersymmetric particles undergo cascade decays to final states with two LSPs 
that give missing energy in the events. 

Our focus is on the signals in neutrino telescopes that should result from the annihila- 



tions of neutralinos that have been gravitationally captured by the Sun 
telescopes are poised to search for the high energy neutrinos of this origin 
experiment at the South Pole is underway and expects to have about 50,000 events from 
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atmospheric neutrinos in the near future 



Sea is currently being built [25, 
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The KM3 detector in the Mediterranean 



30[ |. The IceCube detector has a neutrino 



energy threshold of 50 GeV and the KM3 detector is expected to have sensitivity to neu- 
trinos above a threshold of 10 GeV. These experiments are expected to have the capability 



21| . It is therefore 



to find or limit signals from neutralino DM annihilations in the Sun 
of particular interest to refine the characteristic features imprinted on the neutrino energy 
spectra of this origin. 
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We make substantial improvements on previous theoretical studies 31j by including the 
full spin-dependence of the matrix elements from real and virtual WW, ZZ, Zh, and ti pro- 
duction and subsequent decays to neutrinos, by noting the importance of the spin- dependent 
capture rate in the Sun, and by concentrating on the region of mSUGRA parameter space 
known as the Hyperbolic Branch (HB) 



32 1 or Focus Point (FP) [33| region, where the neu- 
tralino annihilations have the highest rates. (Hereafter we denote this region simply by FP). 
The above mentioned annihilation processes lead to the high energy neutrinos of interest 
for neutrino telescopes. We propagate the neutrinos through the Sun, taking into account 
neutrino oscillations, matter resonance, absorption, and u^- regeneration effects, and give the 
resulting muon-neutrino energy spectra at Earth and the muon energy spectra in km^ area 
detectors. 

Neutralinos are gravitationally trapped and accumulate at the center of the Sun where 
capture and annihilation rates equilibrate over the lifetime of the Sun. The solar capture 



rate of neutralinos in the galactic halo is approximately given by 



C = 3 4 X 1020,-1 Piocai / 270 km/s V / af^ + af, + 0.07< \ / lOOGeV^ ' 
® ■ 0.3 GeV/cm3 V viocai J \ 10~^ pb J \ m^o 

(1) 

where piocai and viocai are the local density and velocity of relic dark matter, respectively. The 
average density is taken to be piocai ~ 0.3 GeV/cm^, but may be enhanced due to caustics 



in the galactic plane 35|. The capture rate is highly dependent on the strength of the 



neutralino interactions with matter 
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361]. The spin-independent (SI) and spin-dependent 
(SD) scattering rates determine how strongly the Sun slows and captures the neutralinos 
and are limited by present direct detection experiments. The factor of 0.07 before a^f comes 
from the relative abundance of helium and hydrogen in the Sun, as well as other dynamical 
and form factor effects For a discussion of the present limits on the SI and SD cross 
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sections from direct detection experiments see Ref. 

As neutralinos accumulate in the solar core, their annihilations deplete the population. 
The competition between neutralino capture and annihilation can be expressed through the 
annihilation rate, F, in the solar core by 

F = ^Cq tanh^ (^^/C^tsun) , (2) 

where Aq = is the annihilation rate times relative velocity per unit volume and tsun = 
4.5 Gy is the age of the Sun. When y/CQAQtsun ^ 1, the capture and annihilation processes 



are in equilibrium, which is expected to be a good approximation for the Sun. Thus the 
annihilation rate is equivalent to half the solar capture rate today since two neutralinos 
create one annihilation event. 

From the thermalization of the captured neutralinos, it has been estimated that the neu- 



tralinos are concentrated with a core region of very small radius, of order O.OlRsun 31| . 



so we can approximate the source of the neutrinos from annihilations as arising from the 



the so^ 
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ar center to the surface, taking 
as well as absorption and 



center of the Sun. We propagate the neutrinos from 
into account neutrino oscillations with MSW effects 
reinjection of neutrinos resulting from neutral current (NC) and charged current (CC) inter- 
actions. Since neutralino annihilation is flavor blind^, the relative populations of different 
neutrino species are the same at the center of the Sun. MSW effects that change the neu- 
trino flavors as they propagate through the Sun are small in the present application because 
the injection spectra are uniformly populated among flavors (i.e. the commutator which 
describes the evolution of the neutrino states as they propagate through the Sun vanishes, 
yielding a steady state solution). In addition, MSW effects are smaller at energies above 
order 100 GeV that are of present interest^. In the results that we present here, we take 
^13 = 0, but we do not anticipate substantial changes for nonzero 6'i3 over its experimentally 



allowed range^. The analysis of Ref. 3l| found that a value of sin^ 26'i3 = 0.04 gave little 
change in the overall neutrino spectra in the high energy region. 

Absorption and reinjection effects due to NC and CC interactions are important. The 
effects from NC interactions are neutrino flavor blind and tend to soften the initial spectra; 
when a neutrino of energy is absorbed, a neutrino of the same flavor but lower energy 
El < E^ is reinjected into the flux. The CC absorption is also flavor blind, but the resulting 
reinjected spectra are not, due to the regenerated Ue and from the leptonic decays of r~ 
leptons. A decays promptly 100% to and 18% to each of z/g and z/^^. Although the Ut 



^ Strictly speaking, secondary neutrinos from the r-leptons and &-quarks from W and Z decays break this 
flavor democracy but these contributions result in softer neutrino energy spectra that are less likely to 
pass the experimental accepjlance cuts. 

^ A concurrent study in Ref. 40| analytically considers the MSW effects of low energy neutrinos from dark 
matter with mass up to 100 GeV. 



Its experimentally allowed range is sin^ 2^13 < 0.19 41 1. 
^ The /^^ leptons can also inject Vf, and z^e, but are stopped by ionizing radiation before they decay, resulting 



in very soft neutrinos 42 1 



4 



number population remains the same, the z/,- energy spectrum is softened by the interactions. 
Additionally, other neutrino flavors receive contributions from the Vr CC interaction. This 
asymmetry in the population of neutrino species introduces small, but nonetheless non- 
negligible vacuum oscillation and NC and CC effects. Once the neutrinos emerge from 
the Sun, they propagate as matter eigenstates to the Earth where they can be detected in 
neutrino telescopes by the high energy muons and electrons produced by the CC interactions 
in the surrounding ice or water. 

The above effects were included in a number of works Q, [ 



3li. In Ref. [31| the 



reconstruction of a general dark matter candidate was studied. However, since the spin 
of the DM particle was not specified, the spin-dependence of the amplitudes could not be 
taken into account. The intent to treat the DM annihilations in a model independent way is 
laudable, but this approach can have serious shortcomings. For instance, DM annihilation 
in the static limit to the two body neutrino-antineutrino process is absent in SUSY models, 
due to helicity suppression. Additionally, processes such as annihilations to weak bosons 
exhibit a striking helicity dependence in the neutrino spectra from their decays. The spectra 
of neutrinos of energy (9(100 GeV) can be substantially altered by the spin dependence of 
the matrix elements. 

In Section [Tll we review the regions of mSUGRA parameter space that can account for 
the observed relic density of neutralinos 43|]. We concentrate on the FP region which solves 



a variety of problems with SUSY phenomenology |44j . l45l . l46l |. In this region of parameter 
space, the lightest neutralino is a bino-higgsino mixture and the neutralino annihilation 
rates are large. In Section UTTl we discuss the neutrino production modes that occur through 
neutralino annihilation in the present epoch. In Section HVl we calculate the neutrino energy 
spectra for each contributing process. We give analytic formulae for the neutrino spectra in 
the static limit where the intermediate state heavy particles are on-shell and show that they 



closely approximate the matrix elements obtained from SMADGRAPH |47l . I48l . I49l . I50l . |51 |. 
Propagation of the high energy neutrinos from the center of the Sun to the km^ area detectors 
on Earth are discussed in Section |Vl The prospects of detecting neutrinos from neutralino 
annihilation are presented in Section IVII We provide our conclusions in Section IVIII 
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II. SUSY PARAMETER SPACE 



In our illustrations we implicitly assume the model of minimal supergravity (mSUGRA). 
In this model the scalar, gaugino and trilinear masses unify at the grand unification (GUT) 
scale. The predictions are determined by the five parameters ^] 



rriQ, mi/2, Aq, tan /5, sign{fi) 



(3) 



Here mo, mi/2 and Aq are the common scalar, gaugino and trilinear masses at the GUT 
scale. The parameter tan (3 is the ratio vacuum expectation values of the up-type and down- 
type neutral Higgs fields and /i is the supersymmetry conserving higgsino mass parameter 
at the electroweak scale. The absence of a Higgs signal from LEP and the relic density 



determination from WMAP 



431 ] limit the mSUGRA parameter space to four distinct regions: 



i) The focus point with large values of rrin and relatively small values of fi can naturally 



yield the observed relic density 
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33 



44 



45 



46 



52 



531]. This region provides 



a decoupling solution to the SUSY flavor and CP problems, suppress proton decay 
rates, and has low fine tuning. The FP has large sfermion masses and mixed higgsino- 
bino dark matter and can be suitable for sparticle mass measurements at the LHC 54] 
. The annihilation modes of the lightest neutralino are predominantly to weak bosons 
and top quark pairs due to the large higgsino-bino mixture of the lightest neutralino. 

(ii) The slepton co- annihilation region has low rriQ and the lightest sle pto n (stau) mass 



is almost degenerate in mass with the lightest neutralino mass 
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The 



resulting co- annihilation can drop the neutralino relic density into the observed range. 



(iii) The bulk region is characterized by low itlq and mi/2 59|, |60|]. In this region, neutralino 
annihilation is dominated by the t-channel exchange of light scalar fermions. However, 
this region is now disfavored by the limits from LEP on the lightest chargino and 



sleptons 
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61]. 



(iv) The Higgs funnel region is characterized by large values of tan/? where 



neutralino annihilates via a light h or broad A resonance 59|, |62, l63|, |64, 



65 



he lig 



itest 
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,l67| 



68]. 



Throughout the remainder of this paper, we adopt the low energy determination of the 



FP region given in Ref. 



461 ] ■ which is preferred for 6 — r Yukawa unification |69|, [70] 



This parameterization adopts tan /5 = 50 and rrit = 174.3 GeV. The central values of this 
narrow wedge are reproduced in Fig. [1] in the space of the parameters fi and Mi at the 
electroweak scale, where Mi is the U{1) gaugino mass parameter. The narrow allowed band 
lies just above the region where radiative electroweak symmetry breaking is not allowed. 
The sfermions and heavy Higgs bosons decouple, so fi and Mi are the critical parameters 
for FP phenomenology of the light gauginos. The mass matrix of the neutralinos in the 
W^, H^) basis given in terms of these parameters is 

/ Ml -giVd/2 giVu/2 \ 

M2 g2Vd/2 -g2Vu/2 

-giVd/2 g2Vd/2 -yu 

\ giVu/2 -g2Vu/2 -/i / 

where gi and g2 are the U{1) and SU{2) gauge couplings, respectively. Gaugino mass 
unification restricts the SU{2) gaugino mass M2 = -^^i ~ 2Mi, and consequently the 
lightest neutralino has substantially more bino content than wino. The mass matrix is 
diagonalized by the rotation matrix Nij 



(4) 



N*M^oN' 



(5) 



The lightest neutralino can then be expressed in the gaugino-higgsino eigenbasis as 



Xi 



° - NnB + N12W + + NuHu. 



(6) 



Using these methods, we numerically diagonalize Eq. (jl]) and determine the mass and 
composition of the lighest neutralino. The composition of the lightest neutralino at some 
representative points in the FP region are given in Table [H 
The mass matrix of the charginos in the {W^, H^) basis is 



M2 \/2Mw cos (3 
V2Mw sin p 



(7) 



The chargino mass matrix is diagonalized by the rotation matrices Uij and Vij 



M2± = U*M^±V^. 



(8) 



The relic density of the lightest neutralino is very sensitive to its gaugino-higgsino content. 
To accommodate the observed relic density, the lightest neutralino is required to have large 



TABLE I: Representative points in the FP region showing the composition of the hghtest neutrahno. 







Ni2 Ni3 


Nu 


90 GeV 


0.80 


-0.20 0.49 


-0.29 


110 GeV 


0.83 


-0.17 0.45 


-0.28 


200 GeV 


0.90 


-0.09 0.34 


-0.24 


400 GeV 


0.87 


-0.06 0.37 


-0.32 


700 GeV 


0.65 


-0.06 0.55 


-0.52 


1000 GeV 


0.24 


-0.04 0.69 


-0.68 




'o 200 400 600 800 1000 1200 1400 200 400 600 800 1000 

Mj (GeV) (GeV) 



FIG. 1: The central FP region in fj, and Mi (left panel) with nit = 174.3 GeV and tan/? = 50 as 
given in [46]. As the mass of the hghtest neutrahno increases, it becomes largely higgsino shown 
by the solid, black curve (right panel). However, the higgsino asymmetry, an important quantity 
for the ZxiXi coupling, decreases as shown by the red, dashed curve. 

bino and higgsino components, leading to the relation /i ^ Mi for light neutralinos as shown 
in Fig. [TJ As the lightest neutrahno mass increases, so do the fi and Mi values and the 
neutrahno mass increasingly controls the annihilation rate until the Xi mass cannot increase 
further as the relic density would increase above the observed value. This behavior is seen by 
the plateau in Fig. [1] at high Mi and yU ~ 1 TeV. The heaviest that the Xi can be in the FP 
region (and still give the correct relic density) corresponds nearly pure higgsino content (i.e. 
^13 + -^14 ~ 1; where N13 and are the higgsino components of the lightest neutrahno). 
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However, since the value of is much larger than the off-block-diagonal elements that are 
0{Mw), the higgsino asymmetry, \Nf^ — Nf^\, that determines the strength of the Zxlxi 
coupling, is very small. 

Since /i ~ Mi throughout most of the FP region, the lightest chargino is dominantly hig- 
gsino and the heavier chargino state is dominantly Wino, with a mass that is approximately 
double that of the lightest chargino. 

III. NEUTRINO SPECTRA 

High energy neutrinos {Ei, > 50 GeV) may be created via neutralino annihilation through 
a variety of channels. Since the FP region has a neutralino with large higgsino and bino 
components, certain diagrams dominate, the most prominent of which are the following: 

(i) Annihilation to weak boson pairs, WW or ZZ, where one or both weak bosons decay 
leptonically. 

Xi W- Xi W- Xi w- 



FIG. 2: Feynman diagrams for the annihilation of the lightest neutralino into WW, ZZ where the 
weak bosons decay into neutrinos. 

The annihilation rate to WW is dominated by the t-channel chargino exchange and 
s-channel Z and Higgs boson exchanges shown in the top panels of Fig. [2j The 
chargino diagram can be enhanced by either large higgsino contents of the chargino 
and neutralino or by a large wino content in the chargino and bino content in the 
neutralino. The ZZ and WW diagrams involving XiXiZ vertex are enhanced by the 




X? 



z 



X? 



z 
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higgsino asymmetry that can be large, as shown in Fig. [TJ For large tan/3, the Higgs 
boson mass is in the decoupling limit where the light Higgs boson mimics the SM Higgs 
boson and the masses of the heavy Higgs states, H^, and if^ are nearly degenerate 
and large. Therefore, the light Higgs exchange processes require both the higgsino and 
bino contents to be substantial in the neutralino exchange diagram. After including 
all contributing processes, it is found that a strong neutrino signal is expected through 
the WW and ZZ channels if M^o > Mw, Mz and < Mi, Ms. 

(ii) Annihilation to Zh, where the Z-boson decays to neutrinos, as shown in Fig. [31 This 
process is enhanced if the higgsino asymmetry is large, due to the XiXi^ coupling 
that occurs in both diagrams. However, a substantial t-channel neutralino exchange 
amplitude requires a large bino component of the lightest neutralino. This process is 
similar to the annihilation to weak bosons except that the emission of the Higgs boson 
changes the helicity structure of the amplitude which alters the shape of the neutrino 
energy spectra. 




X? ^ X? " 

FIG. 3: Feynman diagrams for the annihilation of the lightest neutralino into Zh where the Z-boson 
and decays into neutrinos. 



(iii) Annihilation to top quark pairs as shown in Fig. HI Contributions from annihilations to 
top quarks^, whose primary decay ist^ bW~^, with subsequent W^ — £+z/ decay, are 
dominant if kinematically accessible (M^o > rrit). This process can result in a softer 
neutrino energy spectra compared to direct W^W~ leptonic decays due to the smaller 
available phase space for the neutrinos. The scalar fermion exchange is suppressed 
due to the large rriQ values in the FP region while the Z and Higgs boson exchanges 
require higgsino asymmetry and higgsino-bino dominance, respectively. 



We neglect the QCD corrections [7l| to the annihilation to top quark pairs. 
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h/H/A 




FIG. 4: Feynman diagrams for the annihilation of the Hghtest neutrahno into tt and decays into 
neutrinos. If the scalar fermions decouple, the amplitude for the process in the left panel vanishes. 



(iv) Neutralino annihilation to light quarks and leptons. Due to the Majorana nature of the 
neutralino, the annihilation amplitude to light fermions is suppressed by the fermion 
mass, and direct production of a neutrino- antineutrino pair is therefore negligible^. 
If a 7 or Z boson is emitted, the helicity suppression can be lifted. The leading 
contribution to^Xi ~^ QQ^ via t-channel squark exchange is a^m^o/M| for heavy 



squark masses 72] ■ In addition to squark exchange, s-channel Higgs and Z boson 
processes contribute when the lightest neutralino is dominantly higgsino. We do not 
consider the smaller contributions from these processes. Moreover, since the neutralino 
is above the W mass throughout most of the parameter space in the FP region, these 
modes are negligible^. 

Table [TTl summarizes the dominant subprocess discussed above and gives the correspond- 
ing neutralino annihilation modes and the composition of the Xi that yields a large rate in 
each case. 

The total annihilation cross section must remain relatively fixed in order to provide the 
observed relic density f^^?^^ ~ ^{o-v)^ ' where (erf) is the thermally averaged cross section at 
freeze-out. The DM density is thus inversely proportional to the total annihilation rate to 
neutrinos today^. 



^ Note in the analysis in Ref. [31l | of the neutrino spectra from dark matter annihilation it was assumed 
that the dark matter particle was not Majorana, yielding very clean and distinct neutrino signals from the 
direct dark matter annihilation to neutrino pairs. This possibility is not allowed in SUSY models where 
the neutralino DM candidate is a Majorana particle. 

^ The neutralino lower mass bound is approximately 80 GeV due to the lower mass limit on the lightest 
chargino (M± > 104 GeV) and the gaugino mass unification assumption. 

® Up to factors that include the p-wave terms, which are suppressed in the present epoch where v 10~'^c. 
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TABLE II: Dominant subprocesses for neutrino annihilation to final states that decay to neutrinos. 
The Xi components are given that maximize the subprocess contributions. Note that since Mi ~ fi, 
the Xi 2 3 states have a large higgsino and/or bino fraction and relatively similar masses. The light 
chargino is mainly higgsino while the heavier chargino is mainly Wino. The top quark Yukawa 
coupling is denoted by Yt. The higgsino, bino and Wino contents are denoted by H, B, and W, 
respectively. The H asymmetry is given by {N'^^ — Nf^\ and is shown in Fig. [H 



x?x? ^ WW 


Relevant x? components 


Comments 


t-channel Xi^ 


H 


xf dominantly H 


t-channel xt^ 


W 


X2 dominantly W 


s-channel Z 


H 


enhanced by large H asymmetry 


s-channel h 


H kB 




xlxi ^ zz 




f-channel Xi 


H 


X4 suppressed by small H fraction 


s-channel h 


H k B 




xlxl ^ Zh 




f-channel Xi 


H k {B or W) 


X4 suppressed by small H fraction 


s-channel Z 


H 


enhanced by large H asymmetry 


^ tt 




i-channel t 




Suppressed by M'^o/Mj 


s-channel Z 


H 


enhanced by large H asymmetry 


s-channel h 


H kB 


Enhanced by large Yt 



The annihilation cross section to the neutrinos {a^) via the WW, ZZ, Zh and ti modes 
are shown in Fig. O Once the neutralino mass exceeds the top quark mass, annihilation 
to top-quark pairs dominates, because the Higgs boson coupling to top quarks, Yt, is large. 
However, since the Xi is then dominantly higgsino with large mass, its coupling to the 
lightest Higgs boson decreases since some bino or wino content is necessary. Therefore, at 
TeV scale Xi mass, the annihilation to WW and ZZ dominates. 
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FIG. 5: The annihilation cross section of the lightest neutralino to neutrinos via the WW, ZZ, Zh 
and tt channels in the FP region with tan /? = 50. Note that when the top quark threshold is 
crossed, the cross section for tt production abruptly increases, and the cross sections of all other 
modes correspondingly decrease in order to give the same atot, as is required to maintain the 
observed relic density. 

IV. CALCULATION OF NEUTRINO FLUX 

The high energy neutrino spectra in the FP region are mainly determined by the dia- 
grams in Figs, m [3] and HI We provide analytic results for the neutrino spectra in the cases 
of on-shell WW, ZZ, Zh, and ti and show that these formulas reproduce a Monte Carlo 
(MC) integration of the exact matrix elements provided by SMADGRAPH [49], the super- 
symmetric version of MADGRAPH. In the analytic calculation the phase space integration 
is performed by simple partitioning [t^ and is summarized in Appendix |X1 We calculate 
the spectra for the illustrative points in Table HI 

A. Annihilation to weak bosons 

We calculate the annihilation process XiXi WW* — > Wiu including both on-shell and 
off-shell contributions. The annihilation process of XiXi ~^ W^W~ is isotropic in the static 
limit when the VT-polarizations are not measured. The amplitude of the annihilation to 
on-shell WW is given by the xf exchange with vanishing contributions from the s-channel 
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ZjH diagrams in the static limit, 

M oc e^fi^,e\krre\k<,fk>ikl (9) 

Because the amphtude involves the Levi-Civita antisymmetric tensor, only the two trans- 
versely polarized configurations (++ and ) of W~^W~ are allowed and the longitudinal 

mode is absent. Conservation of CP at tree level implies the two polarization channels carry 
equal weight. 

In the rest frame of W~^, the neutrino u from the decay process ^ i^u is described 
by the normalized angular distribution dN/dcos9 = |(1 — hcos9)^ with h = ±1 for the 
right (left) handed polarized with respect to the axis defined by the momentum in 
the XiXi CM frame. The antineutrino angular distributions can be obtained by a similar 
formula with h ^ —h. 

Including the contributions of both helicities of the W, we find the normalized neutrino 
event rate in the WW CM frame to be 

where (3 = (1 — 4M^/s)^ is the velocity of the W gauge boson and = A/i/4 is the average 
energy. The energy interval A = EdS is such that E^ — A < E,y < Ec + A. The center of 
mass energy of the process is given by s ~ 4M^o in the static limit. The antineutrino from 



W~ has the same energy distribution. The result of Eq. fllOl) is in agreement with Ref. 
Note that the differential rates at the kinematic endpoints of the decay are twice that of 
energy m^o/4. In principle, the distribution in E^ could be useful to roughly determine the 
lightest neutralino mass, provided there is little distortion through NC and CC effects and 
neutrino propagation. Ignoring the W helicity, the differential cross section is instead given 
by a flat E^ spectrum, as often applied incorrectly for neutralino DM. 

In Fig. [6l we show the analytic and numerical evaluation of the neutrino energy spectra 
for XiXi ~^ Wiu in the FP for m^o = 110 and 200 GeV. The analytic result applies only 
above WW threshold. Of note is the clear dependence on the polarization and the slight 
tail on either side of the numerical result, from off-shell decays, in contrast with the sharp 
cutoff of the analytic result. 

The analytic calculation for the neutrino spectra from XiXi ~^ proceeds in a similar 
way, with Z and W couplings and masses appropriately changed, since the neutrinos have 
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FIG. 6: Neutrino energy spectra for XiXi ~^ WW* and WW where at least one W decays lepton- 
ically for m^o = 110 and 200 GeV. Analytic results are given by the dashed curves. 



definite helicity. As the neutralino mass increases, the spectra from WW and ZZ become 
similar since the W/Z mass difference is relatively small. Note that even when one gauge 
boson is off-shell, there is still a rather hard neutrino energy spectra with an appreciable 
rate. This is reminiscent of the significant branching fractions of Higgs boson to gauge 



boson pairs below threshold 75j. Indeed, the neutralino pair in the static limit acts as a 



pseudoscalar state when coupling to fermions 



76| 



B. Annihilation to Zh 

Next we analytically calculate the neutrino energy distribution from the static neutralino 
annihilation, XX ^ ^h, Z vv. The primary process XX ~^ is isotropic in the static 
limit. Only the longitudinal mode of Z is produced in the xX annihilation because of 
conservation of angular momentum. 

In the rest frame of the longitudinal Z, the neutrino v from the decay process Z ^ vv 
is described by the normalized angular distribution with respect to the axis defined by the 
Z momentum in the XiXi CM frame, dN/d cos 6 = |(1 — cos^^). Then the decay of the Z 
boson to neutrinos Z ^ vv gives the neutrino spectra 
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(GeV) (GeV) E^ (GeV) 

FIG. 7: Neutrino energy spectra for XiXi ~^ where the Z decays to neutrinos for m^o = 110, 
200, and 400 GeV. Analytic results are given by the dashed curves. 



where x = SM^oEi,/ (4M^o + ^| — M^) obeys the kinematic range 

\x-l\<p (12) 

We present these analytic results in Fig. [3, along with comparison with our numerical 
results from SMADGRAPH. Since the only longitudinal mode of the Z boson contributes, 
the distribution is very different than that given by the WW or ZZ distribution where only 
the transverse modes contribute. 



C. Annihilation to heavy fermions 

The hardest neutrino energies come from annihilation to top quark pairs, when this 
channel is kinematically accessible. As the threshold for ti is crossed, the cross section to 
ibi^Ui jumps dramatically and dominates the overall neutrino rate. The amplitude of the 
top quark decay, t M^u, has the known structure, {t ■ i^)(b ■ u), which leads to the 
neutrino energy spectra in the rest frame of the t. Boosting into the ti center of mass frame, 
we arrive at the neutrino energy spectra from annihilation to top quarks, assuming a narrow 
VT-width, massless 6-quarks and an average over the t-quark polarization 



1 dN 
NdE^ ~ 

Here we have defined 



(3D 



mm I 



-, — ) I — F f max( 



uo- 



F(y) = (l + 2c.)|-|(|) -0^(1 + 0;) log 



E, 



(13) 



(14) 



where (3 and Et = M^o are the t velocity and energy, respectively, u = m'^/m^, D = 1/6 — 
Lj'^/2+u!^/3, and the expression applies over the energy range (1 — /?) < E^ < {1 + P) -y. 
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FIG. 8: (a) MC result for ofF-shell tt channel, where M o = 170 GeV. (b) Neutrino energy 



distributions for various values of . 



It is interesting to note that the analytic structure of this distribution has two qualitatively 
different neutrino energy distributions, depending on the neutralino mass. Specifically, the 
distribution is divided into three energy regions, the central of which, \Ei, — ^Et{l + u + (3 — 
Puj)\ < jEt\l — uj — P — Puj\, may form a plateau 



1 dN 
NdK 



oc 



Mo 
2 



for M^o > 



2mw 



222 GeV 



(15) 



for M^o < 222 GeV 



These features can clearly be seen in Fig. [Hk, which displays the results for the E^, distribu- 
tions at neutralino masses of 180, 222, and 400 GeV. The high-energy tail of the distribution 
is described by a simple function of the neutrino energy, namely 



dN ( Ey \ ( E^ , \ , 



(16) 



We show MC results for the t*t tht^vi channel in Fig. [8b for the FP with neutralino 



mass m^o = 170 GeV. 

Al 



D. Total neutrino energy spectra 

We combine the above MC results WW, ZZ, Zh, and tt channels to arrive at the neutrino 
energy spectra from neutralino annihilation, as shown in Fig. [H] for neutralino masses m^o = 
90,110,200,400,700 GeV, and 1 TeV. The first two panels illustrate how the WW and 
ZZ spectra change with E^, below the ti threshold. The on-shell ZZ mode turns on at 
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Mz- A characteristic double edge feature is apparent in the total spectra at the 



low and high neutrino energies edges. However, as the neutralino mass is increased, this 
double edge becomes less prominent. Further increasing the neutralino mass includes the 
dominant top spectra. In these cases, the upper edge of the spectra in principle would allow 
a measurement of the lightest neutralino mass. However, as we will see later, propagation 
through the Sun can change the neutrino spectra significantly, especially for large neutrino 
energies, rendering such mass measurements difficult. We verify that the Zh sig nal is 
naturally suppressed; it is considerably smaller than the ZZ and WW modes |36l. |77|. 



V. NEUTRINO PROPAGATION 



The neutrino flux at Earth from annihilating neutralinos in the core of the Sun requires 
a full treatment of the propagation through solar matter and the vacuum between the Sun 
and Earth. To retain the information on coherency of the states, we utilize the quantum 
mechanical evolution equation 



dp 

dr 



-i[H,p\ + 



dp 

dr 



+ 



NC 



dp 

dr 



+ 



cc 



dp 

dr 



4H, [H,p]], 



(17) 



where p is the complex density matrix in the gauge eigenbasis describing the state of the 
neutrino as it propagates. Due to the fact that fast oscillations can only be observed as 
average values due to finite energy resolution, the epsilon term is included in Eq. ( IT7|) 
for convenience in making the energy average jl^ . U^- The choice of e = 0.0 Ir in this 



decoherence term damps all off-diagonal terms to zero 1^. The Hamiltonian, H, includes 
the vacuum oscillation effects from nonzero mass splitting and the MSW terms: 



H 



2E„ 



V2Gi 



N (r) 



(18) 



Here m is the neutrino mass matrix in the gauge eigenstate basis, Ey is the neutrino energy, 
Gp = 1-66 X 10^^ GeV~^ is the Fermi constant and NJr) and Nnir) are the electron and 



neutron densities in the Sun 



7^. 



The source terms, ^ 



NC.CC 



describe the absorption 



and re-injection of neutrinos caused by Neutral Current (NC) and Charged Current (CC) 



processes while the injection source term ^ 
neutralino annihilation in the core of the Sun. 



describes the initial spectra injected by 



«71J 
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FIG. 9: Differential cross-sections dcTu/dEj^ at the solar center for representative neutralino masses 
ranging from 90 GeV to 1000 GeV. Note that when the neutralino mass exceeds the top mass, 
the cross section of the ti channel abruptly increases dramatically, and the other cross sections 
abruptly decrease to maintain the same ^dm- 

A. Vacuum oscillations and MSW effects 



The m)m term in Eq. f|T8|) drives pure vacuum oscillations. In our illustrations we adopt 
the vacuum mass-squared difference values of ^rn^^ = 2.5 x 10~^ eV^, Amf2 = 8.0x 10~^ eV^ 
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and mixing angles 612 = 35.3°, 623 = 45° and 6*13 = 0. The matter interactions with the 
varying solar density suppresses vacuum oscillations. However, since the injected spectra 
are uniformly populated among the neutrino species, vacuum oscillations and MSW effects 
do not change the relative populations, since [H, p] = 0, making the initial state p ~ I3X3 
a steady state. 



B. Source Terms 



The source terms in the propagation equation account for the absorption and re-injection 
of neutrinos as the neutrinos propagate through matter. The CC and NC interactions absorb 
higher energy neutrinos and re-inject them with lower energy. 



1. Injection spectra 



The source term ^ 



dr 



represents the injected initial neutrino spectra from neutralino 



annihilation and is given by ^ 



5{r)5ijj^j^. Here, we make the justified assumption 
that the neutrino flux is initially injected at the center of the Sun 31?]. In the following, 
we illustrate the effect of each term in Eq. f[T7|) . Neutralino annihilation produces neutrinos 
and antineutrinos with each v and v flavor having approximately equal probability. 



2. Neutral Current 



As NC effects are flavor blind, all neutrino species have spectra that are skewed equally 
towards a lower energy. The source term is given by 



dp 

dr 



dE'f^^^i^§^p(K) (19) 

-p{E,) {Np{r)a{uip ^ ulX) + iV„(r)a(z/;n ^ ulX)) . 



The first term represents the re-injection of lowered energy neutrinos with energy E,y from an 
incoming neutrino of energy E^. The negative term accounts for the absorption of neutrinos 
of energy E^. 

In Fig. [TOk we show the neutrino populations in propagation from the center of the Sun, 
where annihilation of a M^o = 200 GeV neutralino occurs, to the solar surface including 
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FIG. 10: (a) Neutrino spectra vs. radius for Ey =50 and 100 GeV including the NC source term, 
(b) Neutrino energy spectra at Earth. Due to the flavor blind NC interactions, the populations of 
the three neutrino flavors remain unchanged as they propagate. 



only the NC effects for = 50 GeV and = 100 GeV. Due to the flavor diagonal NC 
interactions, the three neutrino flavors remain distinct as they propagate. The neutrino 
spectrum of each flavor is equivalent and given in Fig. [TUb . 



3. Charged Current 

The CC effects are also flavor blind to the extent that the r mass is negligible compared 
to E,y. However, the lepton produced in deep inelastic scattering can be stable (in the case 
of the electron), can stop before it decays (muon), or decay before radiation losses decrease 
its energy (tau). The tau lepton is the dominant source of re-injection of neutrinos. The 
channels are 



(20) 



where X denotes hadronic decay products. The source term is 



dp 

dr 



+ 



CC 



dE, 



E: 



^ [Wprr [K) {ET) fr^r [EJE^) (21) 



+U^'^-Prr {ET) fee {ET) ff^e,, {E./ET)], 
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with a similar term for antineutrinos. The = dni^im are flavor projection operators 

where i = 1,2,3 for the e, /x, r leptons, respectively, and F is a 3 x 3 diagonal matrix which 
describes the CC interaction rate for a particular flavor, 

r« = T'cc = (r) cT {uep ^ iX) + iV„ (r) a {u,n ^ ^X) . (22) 

The functions f{u) in Eq. fl2Tl) are the energy distributions of an outgoing neutrino with 
energy Ey, given an incoming neutrino with energy E^J^\ they are calculated in Appendix [D] 
and shown in Fig. [TTl Including this r regeneration effect, the total population does not 
change, but the Vr energy distribution becomes skewed to lower energy. The populations 
^)e,^l energy spectra are increased from the leptonic decays of the r~ lepton, and likewise for 
the for the r+ lepton. 

In Fig. [12] we show the populations of the neutrino similar to those of Fig. [TOl but 
now including only CC interactions. In this case the r decay re-injects r neutrinos and 
occasionally i'e,^i neutrinos, breaking the flavor symmetry. This is clearly seen in the left 
panel of Fig. [T2k for the case of an By = 50 GeV neutrino. 

C. Propagation through the Sun 

We include all the above effects when propagating the neutrino states to the surface of the 
Sun. We illustrate results in Fig. [T3] at neutrino energies of 50 and 100 GeV for neutralino 
masses with values 200 GeV and 1 TeV. Note that as the neutralino mass increases, the 
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FIG. 12: Similar to Fig. 1101 except that only CC interactions are included. The re-injection of z/g, 
from T lepton decays breaks the flavor symmetry among the three neutrino states. 



flavor asymmetry increases. This is due to the larger CC effects at higher energy. 
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FIG. 13: Neutrino spectra after all propagation effects through the Sun are included. The averaging 
of oscillations at larger v/Rq is achieved by the decoherence term in Eq. ()17p . 



D. Propagation to the Earth 

Once propagated to the surface of the Sun, the neutrino states must be propagated in 
their mass eigenstates to the Earth, which is achieved by rotating the density matrix p with 
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the unitary MNSP matrix, V, |80|,l81| 



= vpe,;...v^ (23) 

which diagonahzes the neutrino mass-squared matrix, 

Ml = Vm^mV\ (24) 
where M'^ is the diagonal neutrino mass-squared matrix and m is the neutrino mass matrix 



in the flavor basis. We use the tri-bimaximal mixing matrix for the rotation [8^ . The density 
matrix in the mass eigenbasis is denoted by p. Additionally, the neutrino oscillations are 
averaged by forcing the off-diagonal elements of p to zero in the propagation to Earth. 
This averaging could equally well be achieved by including the decoherence term in the 
propagation equation in Eq. f[T7l) . The neutrino spectra at Earth are quite similar to those 
at the surface of the Sun, as shown in Fig. [1] 



VI. NEUTRINO DETECTION 

Once the neutrinos are propagated to Earth, we need to address their rate of detection in 
a km^ area detector. The fractional rate of annihilation to neutrinos is given by ^— where 
atot includes all SM inclusive modes: WW, ZZ, f f, Zh, hh and a^, = o'xiX°~*'^^ includes all 
the modes discussed above in Section IIIII We calculate the total rates using Calchep and 



can be analytically found in Ref. 



Neutrino Flux at Earth 



We illustrate the detection rates using the simulation for IceCube outlined in Ref. [23|. 
The flux of neutrinos from neutralino annihilation at the Earth is given by 



d^^ _ 1 da^Q^o^^x 1 Cq 



(25) 



dE^ atot dE^ 2 47ri?2' 

where the factor of | is associated with the fact that two neutralinos produce one annihilation 
event; here R = 1.49 x lO^^m is the Earth-Sun distance. The parameter Cq is the neutralino 
solar capture rate of Eq. ([1]) and is dependent on the neutralino mass, velocity and local 
density in the galaxy. 
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FIG. 14: The neutrino energy spectra at production in the Sun, after propogation to the Sun's 
surface, and at the Earth's surface. 



B. FP capture rate predictions 



The solar capture rate of neutralinos in the galactic halo is approximately given by Eqn. 
([1]) and depends on the spin-independent (SI) and spin-dependent (SD) scattering rates^. 

^ The present 90% C.L. limit on the SI cross section from the XENONIO experiment places an upper limit 
on the scattering rate of 8.8 x 10"* pb for a WIMP of mass 100 Gev and 4.5 x 10"*^ pb for a WIMP of 
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FIG. 15: The spin dependent and spin independent scattering cross section of the hghtest neutrahno 
with protons. 



To predict the neutrino flux at Earth, we calculate the spin-dependent cross section for 
the Z exchange^" 

2 



x\-v 

SD 



32/i^J(J 



TT 



,q=u,d,s 



(26) 



M nrri' 



orup 



where J = 1/2 is the spin of the proton, jj,^ = is the reduced mass of the neutrahno 

proton system, T3 is the isospin of quark q and 6g includes the parton distribution function 
and hadronic matrix element of the proton and are given by the values 6u = 0.78, 6d = —0.48 



and 6. = —0.15 
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|. The calculated SD and SI cross sections along the FP region are 
shown in Fig. [THl The SD cross section varies between O{10^^) — O{10^^) pb, while the 
SI cross section is below few x 10^^ pb and is consistent with the limit from XENONIO of 
4.5 X 10^^ pb, making the predicted SD cross section dominate over the typical SI rate by 
a factor of 10^. 



mass 30 GeV 



|. These rates are still higher than the SI rate throughout the FP region, which should 
be probed in the future by CDMS 2007 with an expected reach of 1 x lO^'* pb The SD scattering 
rate off protons, however, is less constrained with a current upper limit of 7 x 10^^ pb from ZEPLIN-II 
for neutrahno masses at 65 GeV [s^. Future detectors such as COUPP can greatly improve these limits 
down to 10~^ pb with a 2 kg chamber of superheated CF3I [87|. However, the SD scattering hmits would 
still be too weak to constrain the FP region. 

Note that the scalar quark exchange is also relevant for the SD scattering rate. However, since the squarks 
have large masses in the FP region, they decouple, and squark exchange contributions are not included 
in the calculation. 
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C. Muon rate at IceCube 



Once the neutrino flux is known, the associated total muon rate through a time T in a 
km^ area detector such as IceCube can be determined by folding the muon production cross 
section with the neutrino flux 



dE„ 



dE„ 



daP{E, 



dE„ 



pp + 



da:{E,^,E,j 



dE„ 



-Pn R^{E^)Aeff (E^) dE,^ + 

(27) 

The densities of protons and neutrons near the detector are taken to be pp = |A^a cm~^ 
and pn = 7:Na cm~^, respectively, where A^^ is Avagadro's number^^. Muons lose energy 



according to 



89| 



dE 

— = -a- PE, 

ax 



(28) 



which can be used to determine the length of a muon track. The parameters a = 2.0 x 10 ^ 



TeV cm^/g and /3 = 4.2 x 10 ^ cm^/g describe the loss rate [22 
range. 



23|, 



89, 



901]. The muon 



R,{E,) = ^\n 



a + pE^ 



a + PEjj^r^' ^^^^ 
is the distance a muon pro pag ates through matter of density p before its energy drops below 
the threshold energy, Ejj^^ 901]. We take Ejj^^ = 50 GeV, which is optimistic for IceCube and 
conservative for KM3. Due to the long muon range, the fiducial volume of the detector can 
be factored into the range and the cross sectional area of the detector, called the effective 
area. The effective area of the detector, A^ff, is calculated for IceCube following Ref. |23| . 
and is given in Fig [161 

The Super-Kamiokande experiment has placed a limit on the flux of muon induced by 
neutrinos from DM annihilations in the Sun of $^ < 5 x 10~^^ cm~^s~^ for a half-angle of 
0-5° 



921 ]. We calculate the induced muon flux according to 



d^ 
dE„ 



dN„, 



(30) 



^eff{Ef,) dEfj^ 

and find a total muon flux above the Super-K bound for 90 and 110 GeV with 1.3 x 
10~^^ cm~^s~^ and 7.9 x 10^^^ cm^^s^^, respectively. The other points we study are well 
below the present bound with 4.5 x 10~^^ cm~^s~^, 2.1 x 10"^'^ cm~^s"^, 3.2 x 10~^^ cm~^s"^ 



Since the muon range is at most 1 km for a 1 TeV muon, the point of muon production can be assumed 
to be in ice, rather than the Earth's crust. 
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FIG. 16: IceCube effective area for given muon energy. 



and 4.4 x 10"^^ cm^^s"^ for 200, 400, 700 and 1000 GeV neutralino masses along the FP 
region, respectively. 

We apply Eq. fl27|) to calculate neutrino signals in IceCube from neutralino annihilations 
as well as backgrounds from atmospheric neutrinos The differential flux of atmospheric 
Ufj^ and Ufi neutrinos, /dE^d cos 6 z, are taken from Ref. jo^, where 6z is the zenith angle. 

To severely reduce the background from atmospheric neutrinos, we include only events 
within a narrow angular cone along the line of sight from the IceCube detector to the Sun. 
We parameterize the zenith angle to of the Sun at the South Pole as 



cos ^7 



cos 



TT 



- + eZf sin (27rr) 



(31) 



where ^2ec^ ~ 23°26' is the maximum declination of the Sun in the celestial sphere throughout 
the year. Here we have denoted the time of year by the dimensionless parameter r which 
spans a full year with the values < r < 1, where r = 0, | correspond to the September 
and March Equinoxes, respectively. The total muon rate along the line of sight to the Sun 
is given by 



dN„, 



dE, 



dE„ 



dr 



M -^E^ JO 

X R^{E,)A,jf{E, 



dal{E, 



Uf,, E^^ 



dEf, 



Pp + 
d cos 9. 



da:{E,^,E^j 



dEy d cos 9 z dr 



dE,. 



R{cos 9z 



-Pn 



(32) 



Backgrounds from solar atmospheric neutrinos also exist, where cosmic rays produce neutrinos in the solar 
atmolphere Q- These backg^unds are expected to be a few per year for I km^ size detector triggering 
on upward going muons. 
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where R{cos6z) = 0.70 — 0.48 cos 6*^ is a detector efficiency for up-going muons that takes 
into account the rock bed below the detector and other angular dependence factors in Ice- 
Cube 231] . The factor of | in the upper limit of the r integration takes into account the time 



that the Sun is below the horizon. The la resolution of IceCube is 1° 95|; we therefore 
include the flux observed along the line of sight to the Sun within a cone of angular diameter 
3°. This dramatically reduces the background from atmospheric neutrinos to that shown in 
Fig.ini 



We show the signal muon event rate and atmospheric background 9J] in IceCube for six 
parameter points in Fig. [T7] and Table lllli Defining the signal region of Ejj^'' < < 300 
GeV, we determine the statistical significance of the six illustrative points in Table lllir ^. 
This choice of the upper cut at 300 GeV is because of the large attenuation of high energy 
neutrinos by the CC and NC interactions in the Sun makes the signal contribution at 
higher energies small. Here, we purposely chose the window to be independent of the 
neutralino mass, but once a signal is seen the window can be narrowed to improve the signal 
to background ratio. In the above window with a few years of data, IceCube can make a 5a 
Xi discovery from Mw up to 400 GeV. The case for 0{Mz) are even better, easily yielding 
a discovery signal in the ffist year of data. 



TABLE III: Muon rates in one year of IceCube data from neutralino annihilations and atmospheric 
neutrino backgrouns and the total statistical significance agtat = S/\/^. The event rates include 
all events where Ejj''' <E^< 300 GeV, with Ef'' = 50 GeV. 
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Note that smearing of the muon energy can yield a tail on the lower end of the spectrum. 
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FIG. 17: Muon event rates and atmospheric backgrounds in IceCube. 



From the signal in Fig. (TTJ one can use the signal rate and the shape of the differential 
muon rate to extract further information on the nature of the neutralino . A determination 
of the neutralino mass can be made from the shape of the muon energy distribution and the 
capture rate of neutralinos in the Sun can be determined from the total signal rate. The 
capture rate in turn gives information on the local galatic DM density and its velocity (cf. 
Eq. ([I])) that can complement inferences from direct detection experiments. 
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D. Limit Spectra 



When the DM particle is very massive and thus produces very energetic neutrinos, the 
features of the neutrino energy distribution tend to be washed out, due to the neutrinos 
interacting more with the solar matter. This gives rise to the so-called "limit spectra," 
which becomes a fair description for E^, greater than 200 GeV or so [31]. The limit spectra 
assumes that (i) oscillations and interactions with matter populate the different neutrino 
flavors uniformly so that the density matrix can be reduced to a scalar, (ii) the re-injection 
spectra from scattering is flat in El/E^, and (iii) the cross section is proportional to the 
energy of the incoming neutrino. In our case, the first assumption is satisfied to zeroth 
order since the initial neutrino populations uniformly populate the three flavors. However, 
r regeneration breaks this uniformity and makes the Ur population different than the 
population^^. The reinjection spectra are preferentially skewed toward lower energies as 
calculated in Appendix [D] and shown in Fig. [11] and are not well approximated by flat 
spectra. The third assumption is valid in the appropriate neutrino energy range of few GeV 
<E^< few TeV. 



For all the limit spectra can be determined analytically [31 1 

p(r,E) = e"^^ |p(r = 0,E) + -y p{r = E')dE' \ (33) 
and approximated as a decaying exponential in the high energy limit 

p(r,E) oc e"^^. (34) 



Here £^ is ~ 100 GeV for neutrinos and ~ 140 GeV for anti-neutrinos [s])]. For very high 
neutrino energies, the spectral distributions become less distinctive and bear little resem- 
blance to the injected spectra. This effect unfortunately limits the ability to detect neutrinos 
of energies greater than 400 GeV or so, as can be seen in FigJTSl Combined with the 
afforementioned low-energy cutoff from the effective area of the detector, we are confined 
to a "window" of possible muon detection energies from roughly 50 GeV to 400 GeV for 
IceCube. 



Oscillations between the Vt and i'^ neutrinos tend to keep these two flavors equally populated. 
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FIG. 18: Given are our exact calculation for neutrinos (solid black line), the approximate spectra 



as described in 



SI], and the limit spectra. 



VII. CONCLUSIONS 

We have made a comprehensive analysis of the prospects for observation in km^ area 
detectors of high energy neutrinos from the annihilations of neutralino dark matter in the 
Sun. The underlying assumptions of this study are the following: 

(i) The lightest stable neutralino of the mSUGRA model accounts for the dark matter 
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density determined by the WMAP experiment 

(ii) The neutrahnos are gravitationally captured by tlie Sun and settle to the core of the 
Sun where they annihilate with the annihilation and capture rates in equilibrium. 
We have assumed the standard estimate of the neutralino capture rate of Eq. ([1]) in 
estimating the signal. 

(iii) The Focus Point /Hyperbolic Branch is the relevant region of the parameter space of 
the mSUGRA model, since the high sfermion masses at the FP solves the SUSY flavor 
and CP problems and gives low fine tuning. 

(iv) The mass of the lightest neutralino has mass above that of the W-boson, so that 
WW, ZZ, and tt are the dominant annihilation channels, giving high energy neutrinos 
from their decays. 

The Focus Point region gives large fluxes of high energy neutrinos for two reasons: 

(i) The spin-dependent cross-section for neutralino capture in the Sun is enhanced by a 
factor of 10^ over the spin-independent cross-section relevant to most current searches 
for direct detection of neutralinos through nuclear recoils. 

(ii) The lightest neutralino is a bino-higgsino admixture which gives a large cross-section 
for the WW and ZZ production processes. 

(iii) The contributions to the neutrino flux are dominated by the top quark contributions 
once the LSP mass exceeds the top quark mass. 

Our calculations included all the physically relevant processes in neutrino propagation 
through the Sun with the quantum mechanical density evolution Eq f|T7j) . including: 

(i) Production at the center of the Sun with equal population of the three neutrino flavors. 

(ii) Oscillations of three neutrinos with vacuum oscillations due to the atmospheric and 
solar mass-squared differences in the tribimaximal approximation to the mixing matrix 
which has a zero ^13 mixing angle and with MSW effects for non-zero ^^13; MSW effects 
are not expected to be significant for high neutrino energy threshold of IceCube. 



33 



(iii) Source terms due to the absorption and re-injection of neutrinos propagating through 
matter, including both CC and NC interactions that re-inject neutrinos with lower 
energy. In particular, we include all neutrino decay products from the production of 
tau-leptons with exact calculations of the tau decay distributions. We find that the 
tau-regeneration effects are not large, but that the absorption effects strongly diminish 
the event rates from neutralinos of mass above a few hundred GeV. 

Our study made major improvements over previous analyses: 

(i) We use take into account the full spin-dependence of the production and decay pro- 
cesses, which have important effects on the neutrino energy distributions. 

(ii) We make quantitative predictions for the absolute neutrino rates, since we specialize 
to the specific FP region. 

(iii) We include off-shell W, Z, and top decays in addition to on-shell. 

(iv) We provide analytic formulas for on-shell decays and show that the numerical cal- 
culations of exact matrix elements obtained with the SMADGRAPH program are in 
agreement. 

We have compared our results with other analyses and approximations where possible: 

(i) We find that the approximation to the density evolution of Eq. f l55]) is reasonably 
good for neutralino masses of order 100 GeV, but shows substantial deviations from 
the exact results at higher neutralino masses: see Fig. [181 

(ii) We find that the limit spectra from Cirelli et al Eq. fl33p is not so good in representing 
the exact results at any neutralino mass, even though it does qualitatively describe 
the fall-off with neutrino energy due to absorption effects. 

(iii) The signal region is broadly defined in muon energy by the 50 GeV threshold cut and 
the rapid decrease at energies above 300 GeV from absorption effects. 

In the detection of muon-neutrino signals from the annihilations, we calculate event rates 
for the IceCube experiment, taking into account 

(i) A muon energy threshold of 50 GeV 
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(ii) The dependence of the detector area on muon energy due to the range of the muons 
in ice. 



(iii) The backgrounds from atmospheric neutrinos as calculated from the flux in Ref . 9J] ; 
we include a 3 degree cone for the angular resolution along the line of sight from the 
Sun to the detector. 

The major conclusions of our study follows: 

(i) The signals from the annihilations of neutralinos in the mass range of 90 GeV to 
400 GeV should be easily identifiable above the atmospheric neutrino backgrounds: 
see Table IIIII About 10 to 200 muons events are created from muon neutrinos by 
neutralino annihilations in the Sun on an atmospheric muon background of about 10 
events (within a 3 degree angular cone the center of the Sun). The breadth of the 
signal is in muon energies a measure of the neutralino mass. The size of the signal is 
governed by the neutralino mass and its couplings to weak bosons and top-quark pairs 
as well as the capture rate of neutralinos in the Sun. 

(ii) The IceCube experiment is most sensitive to neutralinos of mass around 110 GeV, 
above the thresholds for WW and ZZ production by the neutralino annihilation pro- 
cess. However, the signals from neutralinos of mass M^o > rrit, where ti annihilation 
dominates can also be large for neutralinos of mass M^o of order 200 GeV. 

(iii) The absorption effects in the Sun likely preclude the detection in IceCube of neutralino 
annihilations for neutralino masses above 400 GeV. 

(iv) The predicted large signals depend on several assumptions that could be overly op- 
timistic: The Focus Point with large bino-higgsino mixing is the relevant region of 
mSUGRA parameter space (this is subject to validation at the LHC); The Spin- 
Dependent capture cross-section is much larger than the Spin-Independent Capture 
cross-section (this is subject to test in direct detection experiments of WIMP recoils); 
the estimated gravitational capture rate of neutralinos in the Sun is realistic (it could 
be enhanced by caustics or reduced if there is a dark matter density under-density in 
our galactic region). 
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(v) The forthcoming data from IceCube and KM3 offer the prospect for discovery of 
neutrahno DM if nature has cooperated. 
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APPENDIX A: PHASE SPACE PARTITIONING 

The calculations of the cross sections a{xx ~^ WeV) and (j{xX ~^ ZvV) are quite similar. 
Here, we compute the cross section for WW* production, and then give the corresponding 
expressions for the ZZ* process. The cross section (t{xx ~^ WeV) is given by 

1 



where A is the triangular function 



Y,\M^\ds{PSxX^WeV), (Al) 



A(a, 6, c) = + 6^ + - 2{ah + ac + he). (A2) 

The phase space is 

,,iPSxx^WeV)={-) (-j A.(l,^,^)A.(l,^,;;|)*ni. (A3) 
where we integrate over the range mg+m,^ < mx < ^/s—mw- The corresponding calculation 



for ZZ*, with m,, = 0, yields 

2 _2 _2 



d^{PSxx^ZvV)^(^-) j A5(l,-^,-^)dm^, (A4) 
where < mx < ^fs — mz- 



APPENDIX B: NEUTRALINO NUMBER IN THE SOLAR CORE 

It is interesting to estimate the number of the accumulated dark matter particles in the 
solar core. The content of the dark matter located at the core by the equilibrium condition 
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can be modeled by [22 1 



N' = ^ = ^ (Bl) 



Aq {av) 
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The effective volume V of the DM core of the Sun is estimated to be 

= 5.7 X 10^^ cm='(100 GeV/M^o)i (B2) 

The annihilation cross section is dimensionally estimated as 

{av) ^ Ba^/M\ (B3) 

with the coefficient B of order 1. For M^o = 100 GeV, We find ^ 10^^ for the SD capture 
rate. This number is many orders of magnitude larger than the average number N of dark 
matter particles enclosed by a volume in the space of the size of the Sun, iV 4 x 10^°. 

APPENDIX C: NEUTRINO-NUCLEON SCATTERING 

Neutrinos scatter off nucleons via the NC or CC interactions. The cross section expres- 
sions are relatively simple for neutrinos in the energy range we consider {0{1 GeV) < < 
(9(1 TeV)) where the W and Z propagator effects can be neglected since -C M^r. 

1. CC interaction 



For CC scattering, the neutrino-quark cross sections are given by 
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^(^„d lu) = ^{ud^Iu) = (CI) 

dy dy n 

^{vd^lu) = ^(ud^iu) = ^il-yf, (C2) 

dy dy w 

where y = 1 — Ec/E^, = is the CM energy of the subprocess and x is the fraction 

of nucleon momenta imparted on the quark. The resulting scattering off a proton is given 

by H 

^{up^iX) = I dx^{uu ^ id)fu/pix) + ^{i^d iu)fd/p{x) (C3) 
ay Jo "2/ "Z/ 

- ^!!V^:i^(o.i5 + 0.04(1 -y)2), 

7r 

^{up^iX)= f dx^{ud^£u)f^/p{x) + ^{i^u^id)fu/p{x) (C4) 
ay Jq ay ay 

~ ^'""^'^^^ (0.06 + 0.25(1 -1/)^), 

71 
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where fq/N{x) is the parton density function for parton q in nucleon A^. The corresponding 
neutron differential cross sections are [l^ 

^(z/n->£X)= f dx^iuu ^ £d)fu/nix) + ^{ud £u)fd/nix) (C5) 
ay Jq ay ay 

^ 2m^E^^^ + 0.06(1 

TT 

^(z/n IX) = dx'^iyud iu) fa/nix) + ^(^" £d)f u/nix) (C6) 
- ^"'"^'^^^-(0.04 + 0.15(1-2/)^). 



TT 



2. NC interaction 

The NC scattering is calculated in the same way. The subprocess are given by 

^(^g - ^V) = ^(^g - -^'q') = ^ {gl + 9l,ii - yf) , (C7) 
- H') = ^(-g - -'g') = ^ + glii - yf) , (C8) 

where y = 1 — E^/ / Ey and the Z boson couplings are 

12 2 11 1 

gLu = -- -^v^^Ow, gRu = --?>vo?Ow, fi-Ld = -- + - sin^ 6^^/, = - sin^ 6'iy. 

The total differential cross sections for proton target are then given by [1^ 

'^{vp-^v'X) = j\x ^H->^V)/,/p(^) (CIO) 

^ q=u,d,u,d 

~ Tq _^ Q Q22(l - yf] , 

TT 

^{vp^v'X) = j\x J2 ^(^?-^V)/,/p(^) (Cll) 

^ q=u,d,u,d 
71 

and for a neutron target by 

^(un^u'X)= fdx J2 - ^V)/./n(x) (C12) 

g=u,d,u,d 

^ 2m:,E,Gl ^ ^^^^^^ _ ^^2^ ^ 



TT 
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^(z/n ^ u'X) = J dx ^ ^'(l')fg/nix) (C13) 



APPENDIX D: r REGENERATION 

To account for r regeneration, we fold the CC production cross section of leptons with 
decay distribution of the r 

U^^M = nJ'^ (^1 + (<?o. (^) + (^)) ^rf^, (Dl) 

with similar expressions for fp^^^e^ ^^id fu^^p^^- Here, the factors of 1/5 [99!] approximate 
the more precise values in Appendix O In this parameterization, u = E°^^ / E'^J^ and z = 
Et/E'^J^. The normalization factor is chosen such that 

N = / duU^^iu) = { (D3) 
Jo [ 0.18 Ue,f. 

The first factor in each equation arises from the CC production cross section for r~ and r"*", 

respectively. The sum is over each contributing mode (note that for modes with z/g or in 

the final state, there is only one term in the sum), and BEi is the corresponding branching 

fraction. P is the polarization of the decaying r, and is ±1 for r^. g{y) is the energy spectra 

of the neutrinos from the decaying r. We now calculate this spectra for the modes Vr^i^i-, 

VrtjiiVf^l^ and i^rvr, and give a table of all results. 

To calculate the spectra, let /o(x) be the distribution in the static frame of the parent 

r, (e.g. /o(a;) = 2x^(3 — 2x) for i/^ from r). The distribution is transformed into g{y) in the 

fragmentation frame, with y = (x/2)(l + i3cos6), where (3 is the velocity and 6 is the polar 

angle. Each x-bin at Xq or 6{x — Xq) contributes a flat y distribution 9{xo — y)/xo in the 

limit /? ^ 1 for an unpolarized parent. Therefore, 

9o{y) = / dxfo{x)9{x - y)/x = / dxfo{x)/x , (D4) 



Jy 

For Uj. from r, we have 

goiy) = ^ 2x(3 - 2x)dx = ^ - 3^^ + ^y^ (D5) 
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Another way to obtain the result is 

9o{y) = j ^dcos(f) j dxfo{x)5{y - ^x{l + coscj))) 

If the polarization effect in the rest frame is defined by 

1 dN 



(D6) 



dxd cos 6 



fo{x) + fi{x) cos 6 



(D7) 



which, in the laboratory frame, can be parameterized by 89l. Il00 | 

1 dN 



N dy 



9o{y) + Pgiiy), 



then 



9i{y) = + dx{2y - x)fi{x)/x'^ 



For the from r — *• z/^ + /i + t'^ 



= -2x^{2x-l) 



and 



9i{y) 



j 2(l-2x)(2|/-x)rfx= i + ^i/3_3^2 



We now move to the case of in the r decay. The transition probability is 101 1 



Then 



^\M\^ = 64G| (/i- z^r) z^M- {.r-rUrSr) = 16G^m^(l -x)x(l + cos( 
/o(x) = 12x^(1 - x) , fi{x) = +12x^(1 - x) 

g^{y) = [ 12(1 - x)xdx = 2-6y^ + Ay^ 

9i{y) = [ dxl2il-x){2y-x) = -2 + 12y-l8y^ + 8y^ 



(D8) 
(D9) 

(DIO) 
(Dll) 

(D12) 
(D13) 

(D14) 
(D15) 



in agreement with Eq.(96) in 102] • Note that the corresponding spectra for the r"*" decay 
are obtained with the substitution y]^ iy) = ~gl (y). 
For the mode r — >• tt^Ut-, we define 



1 1 
mi — ml 



2mr 



X 



2E,. 



mi 
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TABLE IV: Fragmentation functions for various decay modes of the r lepton, with y = and 

2 / 2 

rx = mj^/m^. 



decay mode 


BF 




^71 (v) 


Z/t-VT 

VrP 
UrX 


0.18 
0.12 
0.13 
0.26 
0.13 


T^^(i -r.-y) 

,\0{1 r,, y) 
T^,0{1 -r,-y) 

03^(0-3 -y) 


i + f y3 - 

2j/-l+rp l-2rp... ^ 
1-rp l+2rp'^'^^ y> 





TABLE V: Fr agme ntation functions for various decay modes of the r lepton, with y 



2 / 2 

rx = mj^/m 



1031]. 



and 



r decay mode 


BF 


50 (y) 


51 (y) 




0.18 


2 - 6y2 + 4y3 


-2 + 12y - 18y2 + gy^ 



The rest frame distribution is a delta function, /o(a;) = S{x — (1 — "Ttt)). We can follow what 
we did before to obtain the boost fragmentation function, 

9o{y) = dxfo{x)/x= dxS{x-{l-r^))/x = e{l-r^-y)/{l-r^) (D16) 



This agrees with formulas by Ref. 103|. For the polarized distribution, the angular distri- 
bution in the parent rest frame is oc 1 — cos 6*, so, 

9iiy) = - [ dx{2y-x)6{x-l+r^)/x^ = -{2y-l + r^)e{l-r^-y)/{l-rl). (D17) 
•^y 
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